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Abstract 

The classical trajectories of single particle motion in a Woods-Saxon and a modified 
Nilsson potential are studied for axial quadrupole deformation. Both cases give rise to 
chaotic behaviour when the deformation in the Woods-Saxon and the 1 2 -term in the 
modified Nilsson potential are turned on. Important similarities, in particular with 
regard to the shortest periodic orbits, have been found. 
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Recent experimental results on metallic clusters reporting abundance variations in mass 
spectra, ionisation potentials, static polarisabilities and collective giant dipole resonances, 
barrier shapes and fragmentation provide us with striking manifestations of shell structure 
effects related to a quantised motion of the valence electrons . The correspondence of the 
electronic shell structure in spherical clusters to the closing of major quantal shells 0, [| [§ 
caused considerable interest in using nuclear shell model type calculations for the descrip- 
tion of metallic clusters || ||, [7[ ||, [|. It turns out that phenomenological potentials used 
traditionally in nuclear physics serve a purpose similar to those obtained within the Kohn- 
Sham density-functional method if the relevant parameters are adjusted appropriately. 
Typical potentials are the Woods-Saxon (and its various modifications) and the modified 
Nilsson potential without spin-orbit term. The considerable lowering of the computational 
time due to their simple analytical form renders an analysis of the stability of large metallic 
clusters feasible. Naturally, the shell numbers have to be larger than the ones used in nuclear 
physics in accordance with the larger number of valence electrons considered for metallic 
clusters. For mesoscopic objects like clusters, deformations, i. e. deviations from spherical 
symmetry of the potential, is as important as in the nuclear physics context. Metallic clus- 
ters can be seen as a 'gift of nature' towards a deeper understanding of the formation of 
shell and supershell structure which is a general feature for any self-consistent theory of 
independent particles moving in an average potential. 

The existence of shell structures is one of the crucial questions addressed by previous 
0, |6], [IT], [7], [| and more recent authors [12|, [13], [L4[|. Obviously, shell structure in the quantum 
mechanical spectrum is associated with periodic orbits in the corresponding classical problem 
IR [TB[. Furthermore, if the corresponding classical problem is nonintegrable and displays 



chaotic behaviour, the shell structure of the corresponding quantum spectrum is affected 
depending on the degree of chaos [T7], |TB|. Since the deformed Woods-Saxon potential as 
well as the modified Nilsson potential are nonintegrable systems, a classical analysis of the 
single particle motion seems to be indicated to shed light on the corresponding quantum 
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mechanical problem. Results of such analysis are presented in this paper. This is relevant 
for two reasons. Since the Woods-Saxon potential and the Nilsson model are both used with 
success in quantum mechanical models for clusters, it is of interest to look at their similarities 
in the classical context; in fact their similarities are not obvious at first glance. Furthermore, 
since the two models show chaotic behaviour as is demonstrated below, it is essential to 
understand whether at least the shortest periodic orbits have similar features; otherwise the 
corresponding quantum problem is unlikely to agree with regard to shell structures. In fact, 
it is the shortest orbits with smallest period that make the most important contribution to 
shell structure in the quantum spectrum fT6||. 

We investigate the classical single particle motion for the Hamilton function (we put the 
mass equal to unity) 

H=l(p 2 r + ^)+Vws(r,0) (1) 

where we use the deformed Woods-Saxon potential Vws( r , 0) = Vo/(l+exp((r— R{6))/d{6))) 
with R = R (l + aP 2 (cos6)) and d = d (l + (VR) 2 /(2R 2 )) pf where P 2 is the second order 
Legendre polynomial and a a deformation parameter. As a second case we consider the 
classical analogue of the quantum mechanical Nilsson Hamiltonian(neglecting the spin-orbit 
term) 

H = \{v\ + V\) + \{^lQ 2 + ^ 2 ) + A(gp, - zp e f. (2) 

In both cases there is cylindrical symmetry. We restrict ourselves to zero value for the 
z-component of the angular momentum, that is p^ = (fir 2 sin 2 9 = 0; we have left out the 
^-dependence altogether. We have chosen cylindrical coordinates g and z in the second case 
while the choice r = y/ g 2 + z' 2 and 6 with tan# = gjz is more appropriate in the first case. 

We first consider qualitatively the case of spherical symmetry which means a = in the 
first and ujj_ = uj z in the second case. Both cases reduce to one degree of freedom since now 
Pe = Or 2 = zp g — gp z is conserved. Closed orbits occur if the radial and angular frequencies 
are commensurate. Rewriting Eq.(2), for u± = u z , as H = l/2(p 2 + pg/r 2 ) + l/2u 2 r 2 + \p 2 e 
we obtain closed orbits if + = n/m with integers n, m. For instance, when n/m = 

1/3, 1/4, . . . the trajectory forms essentially a triangle, a square etc. in the g — z— plane; 
for n/m = 2/5 we get the five star, and so forth. The precise shapes of these geometrical 
figures depend on the magnitude of A in that small values of A produce polygons with 
rounded corners while larger values yield loops at the corners as is illustrated in Fig.(l). 
The appropriate scaling of A is given by the kinematical constraint between the energy and 
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the angular momentum which reads E — Xpg > uj\pe\- Here we emphasise that for negative 
values of A, as used in actual applications, the value of the angular momentum pg is not 
limited in its absolute value for given energy. In particular, the relation implies that (for 
A < 0) if 4|A| > uj 2 / E there is no restriction at all onp^. For the Hamilton function of Eq.(l) 
an appropriate choice of pg can likewise lead to polygon orbits such as a triangle, square, 
pentagon, but also a five star, and so forth. The corners of the polygons are increasingly 
sharp the smaller the value of the diffuseness do (or the larger the value of Rq). The angular 
momentum pg, and therefore the number of corners of the polygons, is now limited by the 
kinematical constraint pg < max r [2r 2 (E — Vws( r 5 0))]. The plain circle (polygon of infinitely 
many corners) is possible only for zero diffuseness where the maximum is reached at r = R . 

Thus a common feature of the two Hamilton functions is, for the trivial spherical case, 
the existence of closed orbits of simple geometrical shapes. Such closed orbits have been 



observed by other authors for the Woods-Saxon potential || [12], [U| but not, to the best of 
our knowledge, for the Nilsson Hamiltonian under consideration. As indicated above, there 
is however a crucial difference between the two: the phase space is non-compact for the 
second case if A < 0. When deformation is invoked both problems become nonintegrable as 
pg is no longer conserved. The symmetrical periodic orbits discussed above are destroyed. 
The onset of chaotic motion can be discerned; for small deformation this happens only in 
parts of phase space. 

For the Woods-Saxon potential we have chosen only a quadrupole deformation of the 
boundary R. The deformation of the diffuseness is a consequence of volume conservation 



IS|. Note that a pure quadrupole deformation of a cavity (do = 0) yields an integrable 



case 



T9I] for the bound state problem. In our case, when a is turned on, the effective 



deformation of Vws(f, 0) implies higher multipoles of order 2 2/ , I = 1, 2, 3, . . . since Vws( r ; Q) 
can be expanded in terms of even order Legendre polynomials. The presence of higher 
multipoles is expected to lead to chaotic behaviour 0] . Results of our analysis confirm this 
expectation. We have solved numerically the canonical equations of motion and obtained 
surfaces of section in the plane g = 0, i. e. the phase space diagrams displayed in Fig. (2) 
are in the z — p z — plane. Our results are general, for demonstration we have chosen the 
parameters Vq = —6 eV, do = 0.74 A, R = 15 A and -3 eV for the energy; this corresponds 
to the Fermi level of a cluster with 300 particles ||. Fig. (2a) represents sections of three 
orbits, the initial values of the one in the centre have been chosen to display the separatrix 
that separates a periodic orbit whose trajectory is displayed in Fig. (3b), and a vibrational 
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mode in the centre of the z —p z ~ plane of Fig. (2a). For this particular value of a the motion 
is still regular in most parts of the accessible phase space. Another periodic orbit with small 
stability islands exists at z = ±17.2 and p z = 0, it is indicated by a pronounced solid dot 
and its trajectory is displayed in Fig. (3a). Only this and the orbit associated with the outer 
part of the separatrix can still be traced when a is increased from 0.1 to 0.16, but only the 
islands surrounding the latter are still significant while the ones at the fringes of phase space 
(z = ±18.7) have virtually shrunk to zero; this is illustrated in Fig. (2b). The onset of chaos 
now occupies larger parts of the phase space. The fast decay of the stability islands for 
only a small range of a- values is remarkable. The centre range in Fig. (2b) is still associated 
with quasiperiodic motion. There are many long time periodic orbits which emerge from 
the centre. The shortest ones are indicated by crosses and their trajectories are illustrated 
in Figs. (3c) and (3d). The most stable and shortest periodic orbit in Fig. (2b) is associated 
with the islands drawn, the trajectory remains that of Fig. (3b). Hard chaos takes over all 
of phase space when a is further increased. The range of a-values considered here is in line 
with values used by other authors |7|] . We mention that omission of the ^-dependence of the 
diffuseness, i. e. putting d = d , does not change the qualitative picture of the surfaces of 
section; it does, however, drastically change an individual chaotic orbit. All orbits displayed 
in Figs. (3) emerge from the centre of the surfaces of section with increasing a. While they 
move outwards in the section with increasing a their corresponding stability region shrinks 
until it disappears. 

For the Nilsson Hamilton function of Eq.(2) deformation is invoked by choosing uj_ > uo z . 
The inequality is associated with prolate deformation. In our context this is the only case 
of importance since oblate deformation is related to prolate by interchanging g and z. 
When A = and u>j_ > uj z the orbits are Lissajou figures, closed orbits are obtained for 
commensurate and uo z . Our interest is directed towards short periodic orbits for A ^ 
where we actually focus on A < because of its physical relevance. We have investigated 
the range of the parameter b, defined by u± = bu z , in the interval 1.5 < b < 2. Fig. (4a) 
displays the phase space structure at g = where the surfaces of section are taken. Note 
that, for A < 0, the whole region to the left and right, i. e. outside the lines z = ±\Jl/ (2|A|), 
is accessible as long as the two lines do not intersect with the ellipse which forms the other 
part of accessible phase space; if |A| is sufficiently large to allow intersection of the lines 
with the ellipse, the phase space becomes connected; then the part of the ellipse which is to 
the left and right of the lines is inaccessible. In Fig. (4b) surfaces of sections for A = —0.01 
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are given for four different orbits. With the choice of energy E = 50 and the frequency 
u z = 7r/2 the two lines are outside the ellipse. Recall that the same pattern is obtained if 
E and A are rescaled such that E • A =constant. With regard to periodic orbits there is a 
remarkable similarity between Figs. (2b) and (4b). In fact, the short periodic orbits are of 
the same geometrical shape and occur in similar regions of the surfaces of section. Compare 
in particular the rather stable orbit of Fig. (4b) with the corresponding orbit in Fig. (2b), 
both have the same trajectory (Fig. (3b)) and their stability islands are situated in a chaotic 
region; further the orbits of Figs. (3c) and (3d) in the centre of the phase diagram which 
forms the stable region of Figs. (2b) and (4b); and the orbit of Fig. (3a) which is situated on 
the fringes of phase space within a tiny region of stability. The diagram of Fig. (4b) refers 
to b = 2, but the pattern as described prevails for 1.5 < b < 2. For larger values of |A| 
hard chaos takes over quickly within the ellipse, in particular when the two lines enter the 
ellipse. However, regular motion prevails outside the ellipse which is the whole area to the 
left of the left line and to the right of the right line. These orbits may attain large values 
of q and z; also the variation of the angular momentum pq is unlimited in principle. In this 
context we stress that the variation of pg ranges typically between —50 and +50 for generic 
orbits inside the ellipse for our choice of parameters. This is significant when considering a 
corresponding quantum mechanical calculation. 

To summarise: the deformed Woods-Saxon potential and the modified Nilsson model 
are both used in quantum mechanical models, more recently in applications to metallic 
clusters. Their classical analogues do not appear to have much in common at first sight, yet 
we have established important similarities between the two. While this provides confidence 
in the quantum mechanical approach where the two models are used for similar purpose, 
our results also call for a certain caution. Both models are potentially chaotic with the 
degree of chaos depending on the deformation parameter in the Woods-Saxon potential 
or on the parameter A in front of the 1 2 -term present in the Nilsson model. Within the 
context of nuclear physics this behaviour is probably of little significance, since only the 
lower end of the spectrum is of interest. For metallic clusters, however, with the much 
larger particle numbers, a higher range of the spectrum becomes relevant; in this region 
chaotic behaviour may well interfere with the search for shell structure. Our results suggest 
that electronic shell structure is not expected to play a major role for the stability and 
formation of metallic clusters if a substantial deformation prevails; however, in previous 



results, where odd order multipoles were considered [18|, the very existence of shell structure 
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was pointed out for prolate deformation. Further investigations have to provide clarity on 
this aspect. Finally we point out that, in view of the large fluctuations in time of the classical 
angular momentum p e , a too severe truncation in the corresponding (stationary) quantum 
mechanical calculation could easily invalidate results - in particular it could conceal the 
onset of chaos in the quantum spectrum. 

RGN gratefully acknowledges financial support from DGICYT of Spain. 
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Figure captions 



Fig.l Typical simple periodic orbits for the spherical Nilsson model. The orbit with the 
loops at the corners is for |A| > uo 2 /(AE) while the other one is for |A| = u 2 /(25E). 

Fig. 2 Surfaces of section for the deformed Woods-Saxon potential for a — 0.1 (top) and 
a — 0.16 (bottom). 

Fig. 3 Typical short periodic orbits for the deformed Woods-Saxon potential in the 
q — z— plane. The same shapes occur in the modified Nilsson model. Figs.(c) and (d) are 
self-tracing orbits. 

Fig. 4 Surfaces of section for the modified Nilsson model. In (a) two situations of acces- 
sible phase space (shaded) are displayed; left: |A| < uj 2 /(AE), right: |A| > uo 2 z /(AE). Four 
orbits for a value of A which corresponds to the left of (a) are displayed in (b). 
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